The Nature of Superfluidity in Ultracold Fermi Gases Near Feshbach Resonances 
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We study the superfluid state of atomic Fermi gases using a BCS-BEC crossover theory. Our approach 
emphasizes non-condensed fermion pairs which strongly hybridize with their (Feshbach-induced) molecular 
boson counterparts. These pairs lead to pseudogap effects above Tc and non-BCS characteristics below. We 
discuss how these effects influence the experimental signatures of superfluidity. 



There has been considerable interest in achieving "reso- 
nance superfluidity" in an ultracold, extremely dilute, trapped 
Fermi gas in which a Feshbach resonance! 1, 2] is used to 
tune the interatomic attraction by variation of a magnetic field. 
These systems allow for the direct study of the crossover prob- 
lem where the nature of the superfluid transition changes from 
BCS to Bose-Einstein condensation (EEC). Recently, several 
groups have observed superfluidity in the EEC regime 
This provides an exciting opportunity for theoretical work 
which can address the entire range of behavior, in anticipa- 
tion of future experiments. 

In this paper we consider the case of a homogeneous gas as 
a first and necessary step. Our goals are to present a theory for 
Tc, for the superfluid phase as well as some signatures of the 
transition, which are shown to be very different from conven- 
tional expectations based on ECS theory. Using current un- 
derstanding of bosonic superfluidity, it is useful to begin with 
a number of inferences about the nature of the superfluidity 
induced when molecular bosons are admixed with fermions. 
For temperatures T above Tc sufficiently strong hybridiza- 
tion (gbl^a^_^a^J between molecular bosons (6^) and fermion 

pairs will result in metastable "pre-formed" pairs (aq_kak)- 
These are necessarily associated with a normal state excita- 
tion gap which represents the energy needed to break them 
apart into fermionic quasi-particles (aj^). For < T < Tc 
the presence of non-condensed molecular bosons similarly in- 
duces the formation of non-condensed fermion pairs. These 
must be present in addition to the usual fermionic single par- 
ticle excitations of the condensate. We may, thus, quite gen- 
erally infer that since the molecular bosons and the fermion 
pairs are so strongly entangled, they must be treated on a sim- 
ilar footing. 

Previous 0, 0] studies of the crossover problem in an 
atomic Fermi gas at T ^ have been based on the work |?] 
of Nozieres and Schmitt-Rink (NSR). This approach, which 
effectively omits non-condensed fermion pairs, addresses Tc 
in a fashion which is not manifestly compatible with the pre- 
sumed ground state [9J. Moreover, pairing correlations are 
only included in the number equation and within an approxi- 
mate and evidently fl^.Tl'l problematic fashion. In this paper 
we address these shortcomings, but similarly base our anal- 
ysis of the standard ECS-like crossover ground state|9], for 
the fermionic component. We begin with the important obser- 



vation that the non-condensed fermion pairs are in chemical 
equilibrium with both the non-condensed molecular bosons 
and the condensate. Consequently, they must satisfy 



T<Tc 



(1) 



This last relation, which is a central equation of this paper, will 
be rewritten shortly in a more concrete form. To implement 
the physical picture discussed above, we follow the theory in 
Ref. 1 12], extending it to include a Feshbach resonance. This 
approach was originally developed to treat high Tc supercon- 
ductors. The physics in the context of boson-fermion models 
is shown here to be different from fermion-only models, al- 
though such hybridi zation models have been also applied to 
the cuprates 1 13, 14]. 

Our Hamiltonian] 1] is given by 
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The sum in a runs over both spin states {T,i}- The free 
dispersion relations for fermions and bosons are given by 
Ek = fcV2m - li and e™ = £^q - 2/i with E'^ = q'^/2M, 
respectively, where we assume /i| = and M = 2m is 
the boson mass. Here f is the detuning of the resonance 
state, C/(k, k') = J/f^SkVk' is the direct s-wave interaction 
and g(k) = gK^k is the Feshbach coupling, with the func- 
tion = exp{ — (fc/i^c)^} providing the momentum cutoff. 
Here we set h — — I- 

The three propagators for fermion pairs, t{Q), molecu- 
lar bosons, D{Q), and single fermions, G{Q) are coupled. 
(Throughout we take the convention J^k = X^k' 
where K, Q, etc. are 4-vectors). With these definitions Eq.Q 
can also be rewritten 



T)-i(O) = =0, T < Tc 



(3) 
H 



The effective pairing interaction is given by 

geffiQ,K,K') = ge//(Q)<^k</'k' with geffiQ) = 
C/ + g2T)o(Q), where T)o(Q) = l/[iQn- E° -ly + 2^I\ is the 



2 



non-interacting molecular boson propagator with Matsubara 
frequency r2„. We have from Eq. Q and the general form of 
the T-matrix, the important result that 



t-i(0)=5e7/(0) + x(0) = 0, 



T<Tr, 



(4) 



We now impose the reasonable constraints that at weak cou- 
pling our results are compatible with the T dependence found 
in the BCS limit for A(T), and for general coupling, the 
fermions are described by the standard T = crossover state 
01 • This, in turn, constrains the fermion pair susceptibility 

xiQ) 



x(Q)-5]G(i^)Go(Q-JO^Lq/2, 

K 

provided also the fermion self energy appearing in G is 



(5) 



(6) 



This form for x ™d E can also be derived by truncating 
the equations of motion for Green's functions at the three- 
particle level. This approach is also closely related to Hartree- 
approximated time dependent Ginzburg-Landau theory 1 15|. 

The T-matrix consists of two contributions: from the con- 
densed (sc) and non-condensed or "pseudogap"-associated 
(pg) pairs. The molecular bosons also contribute to the T- 
matrix via the effective pairing interaction. 



t 



tsciQ) 



9eff{Q) 



i + 9eff{Q)x{Qy 
A? 



T 



-6{Q) 



(7) 

(8) 

(9) 



where A^c = A^c - 9<t>m, with A^c = -C/ X]k(«-kiakT)<y'k 
and (pm — (6q=o)- Without loss of generality, we choose or- 
der parameters Asc and (f)rn to be real and positive with g < 0. 
The order parameter is a linear combination of both paired 
atoms and condensed molecules. These two components are 
connected! 16| by the relation (pm — gl^sc/[{v — 2/i)L/]. This 
implies that A^c = —9eff{0) X]k('^-kiakT)<<5k, as expected. 
Following Eq. (|3}, for T < Tc we may approximate! i3 
Eq. (|6} to yield a BCS-hke self-energy with A^ = A^ 



A2 

pg 



where we define the pseudogap Apg 



(10) 



(11) 



At T = 0, Apg = 0, so that A^c(O) = A(0). Equation 
can be written as 

5e7/(0) + E^^||^^k = 0' T<Tc, (12) 



where £'k = ■\/e^ + A^iyS^. Equation il2l has the form of the 
conventional BCS equation, but the full excitation gap A, as 
distinguished from the order parameter, Asc, appears in the 
dispersion E^. 

We show next that this equation for the excitation gap 
A(T) coincides with the equivalent condition on the molecu- 
lar bosons, given in Eq.(|3|l- Taking the same pair susceptibility 
X in the boson self-energy E^, we obtain 



^BiQ) = -9hiQ)/[i + Uxm 



so that the boson propagator is 



DiQ) 



(13) 



(14) 



After some algebra, Eq.Q leads directly to Eq. il2\ . We may 
now calculate the number equation from the propagators in- 
volved. The number of non-condensed molecular bosons is 
given directly by rif, = — D{Q). For T < Tc, the num- 



ber of fermions is 



rif 



E 



^k ^k 



(15) 



as follows from the condition = 2 J^k ^i-^) "^^^^ a BCS- 
like self energy I^{K). The total number (n) of fermions is 
then 



Uf + 2ni, + 2n" = n , 



(16) 



where ri° 



= (/)^„ is the number of molecular bosons in the 
condensate. 

We now have a closed set of equations for our resonance 
system which require numerical solution. The cutoff function 
(y9k introduces a renormalization of U, g, and v. Extending 
the derivation given in Ref. ! 161 to a nonconstant, separable 
potential, we find that U = TUq, g = Tgo, and i' — i^q + 
ctTg^. Here Uq = ATratg/m < where ai,g is the background 
scattering length and go is the physical scattering parameter 
reflecting the width of the Feshbach resonance. The scaling 
factor r = 1/(1 - aUo) and a = mKrJiin^/^). 

At the physical level, the essential distinction between this 
and previous crossover studies based on the NSR approach 
below Tc is associated with non-condensed pair excita- 
tions lHoil of the superfluid. At a more formal level, it should 
be noted that there is an important difference between Eq. (|5} 
and previous related work !8]. This NSR approach presumes 
that there are two bare Green's functions Gq in xiQ)' the par- 
ticles acquire a self energy from the pairs, but these self energy 
effects are not fed back into the propagator for the pairs. 

To evaluate A(T) n(T) as well as A^dT) and other 
transport properties l! 12l [Tsll below Tc, we solve Eqs. i\ 1> . 
(I12> and il6i for fixed Uq, go, and i/p and a sufficiently large 
cutoff Kc- Due to the divergence of the T-matrix [Eq. @], we 
may Taylor expand the quantity x{Q) in Eqs. il6\ and ( |1 1> to 
first order in f2 and q^: x{Q) = x(0) + ao(«^^« - -Boq^) !15i, 
and expand QeffiQ) and J^b{Q) similarly. This considerably 
simplifies the analysis. 
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FIG. 1: Tc vs detuning vo. Upper inset plots the pseudogap at Tc 
vs fo, and lower inset plots molecular boson contribution to the con- 
densate weight and fermionic chemical potential at T = 0. 
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FIG. 2: Fermionic momentum distribution function at T = and 
T = Tc for the three regimes. PG corresponds to maximal Tc- Inset 
plots the T dependent excitation gap below Tc. 



We begin with calculations of Tc as a function of vq which 
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are plotted in Fig. [l] Here we take go = —4:2Ep/kp and 
Uo = —3Ep/kp as is expected to be reasonably indicative of 
the behavior for Tc in cuiTently trapped atomic Fermi gases. 
For i/Q ~oo only molecular bosons are present and Tc 
approaches the ideal BEC limit. As i^q — » +00, molecular 
bosons become irrelevant and the asymptote of the curve is 
dictated by the behavior of fermions in the presence of U. In 
this paper we have chosen Uq deliberately to be small so that 
the Vq ^ 00 limit is close to BCS. 

Indicated in the upper right inset is the behavior of A{Tc) 
as a function of i^q. The lower left inset shows a plot of 
the molecular Bose condensate weight n^, and the fermionic 
chemical potential /i as a function of i'q. For positive, but 
decreasing i'q, Tc follows the BCS curve until the "pseu- 
dogap" or A{Tc) becomes appreciable. At this point the 
molecular Bose condensate also becomes appreciable. Af- 
ter its maximum (at slightly negative j/q), Tc decreases, as 
does /i, until /i = 0. Beyond this point, towards negative 
i/Q, the system is effectively bosonic. The condensate con- 
sists of two contributions, although the weight of the fermion 
pair component rapidly disappears, as can be inferred from 
the lower inset. Similarly Tc rises, although slowly, towards 
the ideal BEC asymptote, following the inverse effective bo- 
son mass. The corresponding curve based on the NSR ap- 
proach has only one extremum, but nevertheless the overall 
magnitudes are not so different0 01 • There are, however, 
key differences between the behavior of the fermionic exci- 
tation gap A{Tc) and its high Tc counterpart. Because all 
of the condensate comes from molecular bosons in the strict 
BEC limit, A(Tc) reaches a maximum (where the molecular 
bosonic and fermionic weights are comparable), and then de- 
creases, as decreases. This is in contrast to fermion-only 
based models 1 12] where this parameter increases indefinitely 
as the attractive coupling becomes stronger 

In the inset to Fig.|2] we plot the temperature dependence 
of the normalized excitation gap A(T) / A(0) for three values 



of vq/Ef = —200, —5, +200. The second value corresponds 
roughly to the maximum in Tc, where pseudogap effects are 
apparent. The first and last are illustrative of the BEC and 
BCS limits respectively. We, thus, refer below to these three 
values as BEC, PG and BCS cases. While not shown here, 
ji is positive for the latter two and negative for the first. The 
order parameter, Asc, is not plotted here, but for all three cases 
it is rather close to the solid line in the inset. That A(T) is 
relatively constant with T through the superfluid transition is 
to be expected in the presence of "pre-formed pairs", as for 
ua/Ep = -5 and -200. 

It should be stressed that Tc is only apparent in A(T) for 
the BCS case. To underline this point, in the main body of 
Fig. Ewe plot the fermionic momentum distribution function 
rik 1 20], which is the summand in Eq. M5\ . at T = and 
T = Tc. The fact that there is very little change from T = 
to T = Tc makes the important point that this momentum 
distribution function is not a good indicator of phase coherent 
pairing. For the PG case, this, in turn, derives from the fact 
that A(T) is nearly constant. For the BEC limit the excita- 
tion gap, which is dominated by /x, similarly, does not vary 
through Tc- In the BCS regime, A(T) is sufficiently small as 
to be barely perceptible on the scale of the figure. In order 
to address the measurable particle density distribution, these 
observations will have to be incorporated into previous BCS- 
based (i.e., A = A^c) calculations I21L 12211 . albeit generalized 
to include inhomogeneity effects. 

This pseudogap-based phenomenology is weU documented 
in the high Tc superconductors, although for these materi- 
als, penetration depth data (with no analogue here) are di- 
rect probes of the transition to superconductivity. Interest- 
ingly, densities of state measurements in the cuprates also 
show some indications of when order is well established. To 
see how phase coherence enters in the atomic physics con- 
text, we relax the approximation in Eq. ( I10> by noting that 
incoherent or finite momentum pairs (pg) are distinguishable 
from coherent or zero momentum pairs (sc) through life- 
time effects in the self energy. For the pg term we write 



4 




FIG. 3: Fermionic density of states vs energy for the three regimes at 
three indicated temperatures. Note the difference in the scales. 

Epg(w, k) « Apg/{uj + Ck + ij)- For simplicity, we treat 
7 as a phenomenological parameter which is independent of 
T. This distinction between pg and sc is required to arrive at 
general thermodynamic signatures (not discussed here) of the 
transition. 

Figures |3ja)-|3jc) show the resulting density of states 
N{uj) = — 2^j^lmG(ti; + iO,k) and correspond to the 
BEC, PG, and BCS cases, as in Fig. |2l The temperatures 
shown are just above T^, and for T = 0.75Tc and T = O.ST^. 
The BCS case indicates an abrupt transition at Tc, but with a 
very small gap and very low Tc, which may be hard to observe 
experimentally. The BEC case shows very little temperature 
dependence throughout, since the fermions are well gapped 
at all temperatures. Only the PG case, where Tc is maxi- 
mal, indicates the presence of superfluidity, not so much at Tc 
but once superfluid order is well established at T = Q.5Tc, 
through the presence of sharper coherence features, much as 
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seen in the cuprates. 

These plots have important implications for interpreting 
predicted signatures of superfluidity, such as those |23] based 
on laser probing of "atomic Cooper pairs", where it has 
been argued that there is a conceptual analogy between nor- 
mal metal-superconductor tunneling [which measures NiiS)], 
and resonant scattering of laser light. The present work in- 
troduces a caution in interpreting future atomic trap experi- 
ments; Because of the presence of a pseudogap, the signa- 
tures of superfluid onset are not as simple as in BCS or the 
related Bogoliubov-de Gennes theory. In general one has 
to distinguish between the superfluid order parameter and 
the fermionic excitation gap. Nevertheless, superfluid co- 
herence appears to be visible as fine structure effects in the 
fermionic density of states. While the fermionic contributions 
[via N{uj), Uk] do not provide a clear indication of superflu- 
idity, they do establish the nature of the pairing regime: BCS, 
BEC or PG — a "pseudogapped superfluid". 

In addition to investigating and revisiting experimental sig- 
natures, the main theoretical contributions of this paper are to 
establish the presence and role of pre-formed fermion pairs 
at and above Tc, which evolve into non-condensed pairs be- 
low Tc. These pairs hybridize strongly with their molecular 
boson counterparts associated with the Feshbach resonance. 
In this context we present a generalized mean field treatment 
of the broken symmetry phase for ultracold fermionic atoms, 
which unlike other T approaches, connects smoothly to 
the conventional! 9| crossover ground state. Because the Fes- 
hbach resonance has no natural analogue in high Tc systems, 
it remains to be seen whether the differences in these two sys- 
tems will outweigh the similarities. 
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